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Abstract. Numerous ongoing experiments aim at detecting WIMP dark matter 
particles from the galactic halo directly through WIMP-nucleon interactions. 
Once such a detection is established a confirmation of the galactic origin of the 
signal is needed. This requires a direction-sensitive detector. We show that such a 
detector can measure the velocity anisotropy (i of the galactic halo. Cosmological 
N-body simulations predict the dark matter anisotropy to be nonzero, /3 ~ 0.2. 
Baryonic matter has fi = and therefore a detection of a nonzero j3 would be 
strong proof of the fundamental difference between dark and baryonic matter. 
We estimate the sensitivity for various detector configurations using Monte Carlo 
methods and we show that the strongest signal is found in the relatively few high 
recoil energy events. Measuring /3 to the precision of ~ 0.03 will require detecting 
more than 10 4 WIMP events with nuclear recoil energies greater than 100 keV for 
a WIMP mass of 100 GeV and a 32 S target. This number corresponds to ~ 10 6 
events at all energies. We discuss variations with respect to input parameters and 
we show that our method is robust to the presence of backgrounds and discuss 
the possible improved sensitivity for an energy-sensitive detector. 
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1. Introduction 

Several cosmological probes have shown that there is a dark matter component in 
the universe comprising about one quarter of the average energy density [TJ HJ [3]. 
The local density of dark matter in the solar neighborhood is estimated to be in 
the range po = 0.2 — 0.6GeV/cm 3 0]. Several dark matter candidates in the form 
of WIMPs, Weakly Interacting Massive Particles, have been suggested including the 
lightest supersymmetric particle (LSP) in SUSY extensions of the standard model, 
sterile neutrinos and axions - see the review [5J for a discussion. 

Recent years have seen a remarkable progress in the understanding of dark matter 
structures. Universal trends have been identified and quantified using numerical 
cosmological simulations. One of the most discussed general trends lies in the 
behaviour of the universal density profile [6| [TJ |8] , which has likely been explained 
[9] [lOl [IT]. Another general result of cosmological simulations is that the velocity 
anisotropy is zero near the central region, and positive in the outer region |12l 113] . 
There even appears to be a universal connection between the local slope of the density 
profile and the local velocity anisotropy, which allows one to predict that the value 
of the velocity anisotropy near the Earth should be non-zero and of the order 0.2 
[HI HI]. Studies have shown that the velocity anisotropy can have a measurable 
effect on the detection rates of WIMPs [T5J HZ] [TB] . The physical interpretation of 
this velocity anisotropy is that the local dark matter 'temperature' is different in the 
tangential and radial directions with respect to the galactic centre. Thus, a non-zero 
velocity anisotropy of the dark matter presents a sharp contrast with a typical baryonic 
gas. This implies that an eventual measurement of the anisotropy would be a strong 
proof that dark matter really behaves significantly and fundamentally different from 
ordinary matter. This is the main reason why we here present a possible method, 
based on direct detection of dark matter, by which one eventually will be able to 
measure this property characterizing the dark matter halo. 

Detecting dark matter directly involves measuring the recoil of a nucleus which 
is scattered by the WIMP. At present, several direct detection experiments search 
for WIMPs from the galactic halo. Strategies differ, but the best exclusion limits at 
the moment are provided by low background cryogenic detectors detecting phonons 
and either ionization, such as CDMS [H] and EDELWEISS [20 , or scintillation as 
CRESST 21J. However, KIMS [22] are also competitive using only scintillation. A 
promising alternative is dual-phase noble gas detectors measuring scintillation and 
ionization, for example ZEPLIN [33], XENON10 [H] and ArDM [55]. These are 
relatively easily scalable to ton-mass detectors. The only collaboration to claim a 
detection so far is DAMA [53] which relies on detecting the weak annual modulation 
of the signal rate induced by the motion of the Earth [37]. However, the claim appears 
to be ruled out by CDMS [H] and is heavily disputed. Of particular interest to the 
present work are the direction sensitive detectors DRIFT [28] and NEWAGE [2"5] . 
DRIFT is a 1 m 3 negative ion time projection chamber (TPC) situated in the Boulby 
Mine in the UK. The collaboration has provided proof-of-principle and are running 
the second stage of the detector. NEWAGE is a micro-TPC in the R&D stage. 

Direction-sensitive detectors search for a WIMP signal induced by the solar 
motion through the WIMP halo [30 . This causes a large forward-backward asymmetry 
in the recoil signal rate as shown in figure [T] The magnitude of the asymmetry depends 
mainly on the solar orbital speed. Analyses suggest that less than 10 WIMP events 
in a direction-sensitive detector may confirm the signal as being galactic due to the 



The dark matter halo velocity anisotropy 



3 




Figure 1. Hammer- Aitoff projections of the sky in galactic coordinates showing 
the directional nuclear recoil rate in our fiducial model. On the left all recoils 
are shown, on the right only those with a recoil energy E ICC > 100 kcV. The Sun 
moves towards 90° longitude on the Galactic equator. 

large forward-backward asymmetry [3TJ [35J |33J |M] • Hence, while direction sensitive 
detectors are less likely to provide first detection, they may well provide a crucial 
confirmation of the galactic origin of a WIMP signal. However, the asymmetry is also 
weakly dependent on the velocity anisotropy of the dark matter halo and therefore a 
careful measurement of the asymmetry allows for a determination of f3. 



2. Modelling 

In this section we discuss the components of the Monte Carlo simulation. The 
simulation generates a detector event by randomly selecting a WIMP velocity in the 
galactic halo, transforming the velocity to the detector and calculating the detector 
response to the event. 



2.1. Dark matter velocity distribution 

We make the simplest assumption for the velocity distribution which is a modified 
Maxwell-Boltzmann distribution , 

in which the one-dimensional radial and tangential velocity dispersions are related 
through (3 = 1 — o~1 ja^.. As explained in the introduction, numerical simulations 
predict (3 ~ 0.2 near the Earth. When varying (3 we fix the total velocity dispersion 
a = \/of + 2of at a — v — 230 kms -1 , the circular orbit velocity. The distribution is 
cut off at an escape velocity of v csc = 600 kms -1 . More general distributions [35J 
will require a more detailed study. 

The velocity of the Earth through the Galaxy is calculated as specified in 
Appendix B of [37] to obtain the WIMP velocity v in the detector. The velocity 
is the sum of the Earth's orbit around the Sun ve, the solar motion vq with respect 
to the local standard of rest and the galactic orbital speed of the local standard of 
rest, which is just the circular orbit velocity vq. The solar motion measured by the 
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Hipparcos satellite [38] is, in Galactic coordinates, v Q — (10.0, 5.2, 7.2) kms 1 . We 
refer to [37 for the calculation of ve as a function of time and date. 

2.2. WIMP-nucleus interaction 

A WIMP can be detected by the nuclear recoil produced when it scatters off a target 
nucleus in a detector. We consider the case where the dominant channel is spin- 
independent clastic scattering, which is coherently enhanced by the number of nucleons 
at low energies. The cross section depends on the momentum transfer q as 



where F(q) is the nuclear form factor and ao is the cross section in the limit of zero 
momentum transfer. Following [37], the WIMP-nucleus interaction is modelled by the 
Helm form factor [55] . 



thickness parameter. For elastic scattering the momentum transfer is given by 



where /i — m t mw / '(rnt+rnw) is the reduced mass of the WIMP-target system, v is the 
laboratory speed of the incoming WIMP and 9 is the recoil angle of the nucleus with 
respect to v. For each event, the maximum possible momentum transfer is determined 
by the velocity of the WIMP. The distribution |F(g)| 2 is then randomly sampled in the 
interval up to the maximum momentum transfer to determine the actual momentum 
transfer. This fixes the scattering angle 9 through Q. 

2.3. Detector 

The directional WIMP signal is unavoidably smeared out by the fact that only the 
nuclear recoil is observed, not the WIMP itself. In addition to this, the detector will 
have a limited angular resolution with which the initial nuclear recoil direction can be 
reconstructed. For example, in the DRIFT time projection chamber, the ionization 
cloud from a typical recoil will drift onto only a few anode wires in the readout [40 . 
This, coupled with the charge diffusion, necessarily limits the accuracy with which 
the recoil can be reconstructed. High energy recoil events produce longer ionization 
tracks, hence it is easier to measure the direction of more energetic events. 

In an actual experiment the angular resolution must be carefully measured. Here, 
we model it by rotating the recoil velocity in a random direction by an angle a drawn 
from the Fisher distribution on the sphere p(a) oc exp(«;cosa) [41] . The parameter n 
fixes the width of the distribution with larger k corresponding to a more centralized 
distribution. We consider perfect reconstruction of the direction, k = oo, as well as 
k — 5 and n — 2.3 corresponding to half the sampled angles being greater than 30° 
and 45°, respectively. 

Another parameter characterizing the detector is the detector threshold energy 
Ethr- A realistic detector cannot measure the direction of nuclear recoils below this 
energy since low energy recoils will not move sufficiently long distances in the detector. 
Naturally the total signal rate is lowered for increasing threshold energy. 



a(q) = a \F(q)\ 2 , 



(2) 




q = 2[j,vcos9, 



(4) 
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Figure 2. Distributions of nuclear recoils with respect to galactic longitude, 
normalized to unity, for three values of /3. An energy threshold of Bthr = 100 keV 
has been applied, other parameters are as in the fiducial model. 



2.4- Calibrating the asymmetry-anisotropy relation 

The observed asymmetry depends on (3 in a non-trivial way which necessitates Monte 
Carlo calibration. 

We define the angle 4> between the observed nuclear recoil direction and the 
direction of the solar motion through the halo. If (f> is less than a chosen acceptance 
angle a acc the event is counted as a forward event while if cj> is greater than 180° — a acc 
the event is a backward event. The observed asymmetry is the difference between 
the number of forward and backward events. From figure [2] it is evident that the 
asymmetry is large, compared to the total number of signals, but the dependence on 
f3 is weak. The relative asymmetry is the difference between the number of forward 
and backward events with E rcc > E t hi, divided by the total number of generated 
events, regardless of recoil energy. In other words, it is the probability that a random 
detector event will add to the asymmetry given E t hr and a acc . 

The relation between the relative asymmetry and the velocity anisotropy is 
calibrated by simulating 10 8 events for (3 in the range (—0.5,0.5) in steps of 0.05. 
For each f3, the relative asymmetry is tabulated for a number of threshold energies 
and acceptance angles. Third order polynomials are fitted to the calculated relative 
asymmetries and used as calibration curves. Figure [3] shows the calibration curves 
at threshold energies of OkeV and 100 keV. Note that for low threshold energy, the 
relative asymmetry increases with (3 while the opposite is true for high threshold 
energy. This change in behaviour occurs in the region of -Ethr — 40keV, increasing 
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Figure 3. Calibration curves for (a) E t h r = OkeV, (b) E t h r = 20keV, (b) 
^thr = 40keV and (d) -E t i lr = 100 keV, acceptance angles as labelled. The 
calibration is for the fiducial model, i.e. k = oo, my/ = 100 GeV and 32 S target. 



with the smaller acceptance angles. 

The calibration depends on both the WIMP and target masses, mw and m t , 
on the assumed velocity dispersion a, and on the angular resolution of the detector, 
k. Further, each configuration of detection energy threshold and acceptance angle 
results in a different calibration. We calculate calibration curves for £7 t hr/(keV) G 
[0, 20, 40, 70, 100, 140, 180] and a a cc € [90°, 60°, 45°, 36°, 30°, 26°, 22°]. 

3. Results 

Now we discuss the sensitivity, i.e. the mean or expected accuracy, with which the 
velocity anisotropy can be measured, depending on the experimental configuration. 
We simulate a large number of experiments, each measuring from a number of 
observed nuclear recoils. The measurement is done by converting the observed relative 
asymmetry in each experiment to a value for (3 using the calibration curves discussed 
above. This yields a distribution of measurements and the sensitivity is the width of 
this distribution. 

3.1. Signal rates 

We do not calculate specific rates, rather we assume a WIMP mass and a target 
nucleus and estimate the sensitivity from a number of generated events for different 
i?thr and a acc . Increasing E thr lowers the signal rate by a factor depending on the 
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Table 1. Fraction of nuclear recoils with recoil energy above -E t hr- 



E thl (keV) 







20 


40 


70 


100 


140 


180 


32 S target, raw 


= 100 GeV 


1.0 


0.35 


0.18 


0.074 


0.030 


0.0082 


0.0018 


12 C target, raw 


= 100 GeV 


1.0 


0.22 


0.073 


0.013 


0.0014 






32 S target, mw 


= 500 GeV 


1.0 


0.44 


0.27 


0.14 


0.073 


0.031 


0.013 



masses. Table [T] lists the fraction of recoils with energy above a given threshold for 
some combinations of mw and target nucleus. For example, out of 1000 recoils only 
about 30 would have recoil energy greater than 100 keV if mw = 100 GeV and the 
target is 32 S. It is therefore natural to compare sensitivities for different threshold 
energies between simulations with the same number of total generated events. 

The actual signal rate in a detector is proportional to po/mw, where po is the 
WIMP density in the vicinity of the solar system. 

3.2. Fiducial model 

We consider a fiducial model with velocity anisotropy (3 — 0.1, a WIMP mass of 
mw = 100 GeV, a 32 S target nucleus and perfect reconstruction of the recoil direction, 
k = oo . We simulate a total of 10 9 detector events which we group as 1000 experiments 
obtaining 10 6 WIMP detections each. The distribution of measured /3's for these 1000 
simulated experiments is shown in figure [4] for E^ — 100 keV and a acc = 45° . The 
width up of the distribution is the desired estimate of the experimental sensitivity to 
(3. Explicitly, a detector reconstructing the direction of the roughly 3 x 10 4 events with 
recoil energy greater than 100 keV (out of the 10 6 events at all energies) is expected 
to be able to measure (3 with an accuracy of ap = 0.029, if (3 = 0.1. 

The sensitivity in the fiducial model is shown in figure [5] for various acceptance 
angles as a function of threshold energy. The best sensitivity is obtained for the above- 
mentioned configuration, a acc = 45° and E t \ u - = 100 keV, for which ap = 0.029. Very 
similar results are obtained for £thr = keV while intermediate threshold energies are 
less optimal. The dependence on a acc is weak for large threshold energies, but rather 
strong at low and intermediate E^ T . The optimal threshold energy at 100 keV should 
be understood as the best compromise between a steep calibration curve and a large 
number of events. Explicitly, the asymmetry is more enhanced if only high energy 
events are detected, but the lower number of events reduces the signal to noise ratio. 
For zero threshold energy, the opposite happens as there is a large number of events 
but each event contains less information. The poor performance at intermediate £thr 
is mainly due to the calibration curves changing from positive to negative slope as the 
threshold energy increases while the increase in ap for E t hr > 100 keV is due to the low 
number of events with sufficient recoil energy. The simulated experiments reproduce 
the input (3 = 0.1 consistently. The only unfortunate exceptions are for a acc = 90° 
and i?thr = 20 or 40 keV in which case the flatness of the calibration curves and the 
low signal per event smear the distribution out over a wide range of /3's. 

The measured value of the relative asymmetry in the 1000 experiments is close 
to Gaussian since the number of forward events is very large, the number of backward 
events is very small and each number is Poisson distributed. However, since the 
calibration curves are not straight lines the distribution of measured /3's is distorted. 
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Figure 4. Binned distribution of measured [3's in the fiducial model with 10 
detector events, Ethr = fOOkeV and aacc = 45°. The width of the distribution 
defines the sensitivity for this experimental configuration. The curve shows the 
Gaussian fit which has = 0.029. 



7=-3 = Tv^ ( 5 ) 



The measure of this distortion is the sample skewness, which is the ratio of the third 
sample moment to the second, 

[£,(&- OT/ 2 ' 

In particular, we find large, negative skewness 7^—1 for low threshold energies as 
well as for a acc = 90°. In these cases a Gaussian fit is not a good representation of 
the distribution of /3's and we take the sample standard deviation of the unbinned 
distribution as the sensitivity instead. 

The dependence on the number of detected events is investigated by regrouping 
the 10 9 simulated events into more experiments with fewer detections. Figure [6] 
shows the obtained accuracy for 4 x 10 5 events, i.e. 2.5 times fewer than the fiducial 
configuration. The best 073 is again obtained at -Ethr = 100 keV and a acc = 45°, for 
which 1J/3 = 0.047. In general we find that ap cx 1/yN for experiments with at least 
10 5 events, as expected. However, the skewness becomes more and more pronounced 
as the number of detector events is reduced. 

In summary, to make a coarse measurement of the velocity anisotropy of the dark 
matter halo with an accuracy of about 0.1, about 10 5 WIMPs at all energies would be 
needed. However, a detector which is insensitive to the recoil energy should ideally only 
count the nuclear recoils with energy greater than 100 keV. A precision measurement 
would require more than 10 6 events, corresponding to a ton-scale direction-sensitive 




Figure 5. Sensitivity crp of measurement of [3 as a function of detector energy 
threshold for various acceptance angles in the fiducial model. 



detector operating for several years if the cross-section is 10 _7 pb. These numbers, 
compared to the present status of WIMP directional detection experiments, are 
of course very large and would require a large dedicated experimental programme. 
Additionally, a cross-section of 10 -7 pb may be too optimistic - for example, a recent 
study found that the most favored region would be the range 1CP 10 — 10 _8 pb [42]. 

3.3. Parameter influence 

Next, we investigate the consequences of varying each of the parameters of the fiducial 
model independently. 

• First and most important, a limited detector angular resolution, k = 5, is found 
to increase ap by 30%. A poorer detector with n — 2.3 increases ap by 75%. 
The optimal acceptance angle in both cases is a acc = 60°. In these cases, the 
073's obtained at zero threshold energy are notably poorer than at 100 keV, unlike 
in the fiducial model. One might expect the resolution would have a stronger 
effect but it should be remembered that the WIMP recoil distribution is already 
smeared out by the nuclear recoil distribution so the relative decrease due to a 
finite k is small. 

• If the target nucleus mass is lowered, the reduced mass of the WIMP-nucleus 
system is decreased and, from Q, so is the width of the recoil angle distribution. 
Hence the nuclear recoil direction resembles the incoming WIMP direction better. 
For a 12 C target, 073 is decreased by 15%. The lower reduced mass also lowers 
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Figure 6. Same as figure[5]but with 4 X 10 5 events per experiment. 

the optimal threshold energy to E t hr = 70 keV. 

• If the WIMP mass is increased, the reduced mass is also increased and the 
recoil angle distribution becomes wider, resulting in poorer sensitivity. For 
mw — 500 GeV we find ap is increased by 10% and the optimal threshold 
energy is 180 keV. It should be noted here that since the WIMP flux is inversely 
proportional to mw, it would take five times as long to record the same number 
of events if mw — 500 GeV rather than 100 GeV. 

• Finally, if the 'true' value of (3 is instead —0.1 the smaller slope of the calibration 
curve causes ap to increase by about 50%. The opposite is true if (3 = 0.3, in which 
case (7/3 decreases by 40%. Hence, if the actual value of (3 is larger than the 0.1 
in our reference model, the number of events needed to measure (3 is significantly 
smaller. For example, measuring (3 = 0.3 to a precision of about ±0.06 would 
require roughly 3 000 events above 100 keV. Figure [7] shows the variation of the 
sensitivity with the assumed true value of (3. 

The best obtained ap for these parameters and the corresponding optimal acceptance 
angle and threshold energy are summarized in table [2j 

3.4- Background 

So far we have assumed a zero background level. Now we discuss the influence of 
a nonzero background which is assumed to be isotropic in the Galactic frame when 
averaged over time [43] . We consider the impact on the optimal detector configuration 
for the fiducial model, £W = 100 keV and a acc — 45°. For each simulated experiment 




Figure 7. Variation of the sensitivty with the input value of j3 for a a cc = 45° 
and 60° and -Ethr = 100 keV. The dependence is primarily determined by the 
calibration curve. 



Table 2. Impact of varying the parameters of the fiducial model independently. 



Parameter 


Best ap 


E th r (keV) 


^acc 


Fiducial model 


0.029 


100 


45° 


K = 5 


0.037 


100 


60° 


k = 2.3 


0.051 


100 


60° 


m w = 500 GeV 


0.033 


180 


45° 


12 C target 


0.025 


70 


45° 


= -0.1 


0.045 


100 


45° 



we add a number of background events to the forward and backward signal events. 
These two numbers are drawn from a Poisson distribution with mean equal to a 
fraction of the number of forward signal events. The asymmetry is then calculated 
as before and the calibration curve is used to determine the measured value of (3. 
We find that the degradation due to nonzero background is benign, as long as the 
signal is not weaker than the background. For a 25% background level the sensitivity 
<jp = 0.037, an increase of less than 30%. For a 100% background, <jp = 0.06 while 
for a 400% background 073 = 0.10. The distribution of measured /3's still reproduce 
(3 = 0.1 consistently as shown in figure [8] 

The well-behaved performance of the sensitivity with respect to nonzero 
backgrounds can be attributed to our method of measuring the velocity anisotropy 




Figure 8. Binned distributions of measured [3's in the fiducial model as in figure 
[4] but including a background level set relative to the number of forward events. 
The sensitivities are = 0.037 and an = 0.06 for 25% and 100% background 
levels, respectively. 



from the difference of the number of forward and backward signals, defined within 
equal solid angles. This method is largely insensitive to the background level, since 
the expected number of background events is the same in the forward and backward 
bins. 



3.5. Energy resolution 

We have not assumed that the detector can measure the nuclear recoil energy. 
However, as is evident from the calibration curves in figure [3] the forward-backward 
asymmetry increases with (3 for low energy recoils but decreases for high energy recoils. 
Hence a 'cleaner' signal can be obtained if the experiment is able to place events into 
energy bins with some resolution. 

We have investigated the possible improvements by binning events into 15keV 
bins according to their recoil energy, i.e. with no additional detector resolution effects. 
For each energy bin and acceptance angle we calculate separate calibration curves. 
The flip from positive to negative slope of the calibration curves takes place at about 
-E-rcc = 60kcV for the fiducial model. Hence, at intermediate recoil energies there 
is virtually no sensitivity to f3. Following the usual procedure, we calculate the 
distribution of measured /3's in each energy bin. For 10 6 events distributed over all 
bins, the measurements reproduce the input (3 for low (< 50 keV) and high (> 100 keV) 
recoil energies. The best sensitivity in individual bins is obtained in the lowest energy 
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bin, E rcc E (0, 15) keV for which 073, .e = 0.032. If the results of the bins that reproduce 
= 0.1 are combined through a^ 2 = J2i a ~a a sensitivity of cxg = 0.20 is achieved. 
This corresponds roughly to the improvement obtained by our standard detector with 
no energy resolution taking data for twice as long. 

It is interesting to note that a direction-sensitive detector with high energy 
resolution will be able to extend the study of [44] to extract the shape of the full 
3-dimensional velocity distribution function. 

4. Discussion 

We have investigated the possibility of measuring the velocity anisotropy (3 of the 
galactic dark matter halo in a direction sensitive WIMP detector using Monte Carlo 
methods. The measurement is based on the well known forward-backward asymmetry 
in the directional spectrum of WIMP-induced nuclear recoils. A non-zero f3 alters the 
magnitude of the asymmetry slightly which makes it possible to measure f3. 

We find that in excess of 10 5 events across all energies are needed to make a coarse 
measurement of (3. An experiment measuring 3 x 10 4 events with recoil energy greater 
than 100 keV, equivalent to 10 6 events at all energies, should be able to measure f3 to 
a precision of 0.03. This result is obtained for an acceptance angle of a acc — 45°, a 
velocity dispersion of a = 230 kms -1 , a WIMP mass of 100 GeV and a 32 S target. Such 
a measurement would provide a strong proof that dark matter behaves fundamentally 
different from baryonic matter. We also note that if j3 is actually larger than the 0.1 
we have assumed, the required number of events can be reduced significantly. 

We have investigated the dependence of the sensitivity with respect to the angular 
resolution and the threshold energy of the detector, the masses of the WIMP and the 
target and the background level. The main point is that low energy and high energy 
nuclear recoils carry more information about the velocity anisotropy, implying a steep 
calibration curve, while intermediate energy recoils show little sensitivity. Since the 
calibration curves for low and high energy recoils have positive and negative slope, 
respectively, an energy-sensitive detector may improve the sensitivity. 
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